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ABSTRACT 

Evidence is accumulating suggesting that collisionlcss shocks in supernova remnants (SNRs) can 
amplify the interstellar magnetic field to hundreds of microgauss or even milli-gauss levels, as recently 
claimed for SNR RX J1713. 7-3946. If these fields exist, they are almost certainly created by magnetic 
field amplification (MFA) associated with the efficient production of cosmic rays by diffusive shock 
acceleration (DSA) and their existence strengthens the case for SNRs being the primary source of 
galactic cosmic ray ions to the 'knee' and beyond. However, the high magnetic field values in SNRs 
are obtained exclusively from the interpretation of observations of radiation from relativistic electrons 
and if MFA via nonlinear DSA produces these fields the magnetic field that determines the maximum 
ion energy will be substantially less than the field that determines the maximum electron energy. 
We use results of a steady-state Monte Carlo simulation to show how nonlinear effects from efficient 
cosmic ray production and MFA reduce the maximum energy of protons relative to what would be 
expected from test-particle acceleration. 

Subject headings: acceleration of particles, shock waves, cosmic rays, supernova remnants, magnetic 
field, maximal energy 



1. INTRODUCTION 

Magnetic fields with strengths well above what can 
be expected from shock compression of the interstel- 
lar medium (ISM) field, £?ism, are inferred in super- 
nova remnants (SNRs) from four types of observations: 
(1) sharp X-ray edges or filaments (e.g.. iBamba et all 
12001 iWarren et all l2005t iCassam-Chenai et all 12007^ 



(2) broadband fits of emission from radio to 7-ray 
energies (e.g., [Cow sik & Sarkar I98(| JVink fc Laming! 
l2003 HVolk. Berezhko fc Ksenofontovll2005l). f3) sp ectral 
curvature in radio emission (|Revnolds fc Ellisonl [1992T) 
and, (4) with the recently reported Chandra obser- 
yations of brigh t filam ents in SNR RX JI7I3. 7-3946 
(|Uchivama et al.l l2007f) . rapid time variations of non- 
thermal X-ray emission. The most likely explanation 
for these fields is that -Bism is amplified as part of 
the efficient production of cosmic rays by nonlinear 
diffusive shock acceleration (DSA) (e.g.. iBell fc Lucek 
200lMAmato fc Blasill2dollVladimirov. Ellison fc Bvkov 
20061 ). Since the magnetic field largely determines the 



maximum particle energy a shock can produce, magnetic 
field amplification (MFA) is critical for interpreting GeV- 
TeV emission and determining if SNRs produce cosmic 
rays to the 'knee' and above. Equally important is the 
fact that if strong MFA is an intrinsic part of nonlin- 
ear DSA, unexpectedly strong magnetic fields may be 
present in other systems undergoing efficient shock accel- 
eration such as radio jets and gamma-ray burst (GRB) 
afterglows. The detailed modeling of broadband observa- 
tions of SNRs is the only practical way to determine the 
importance of nonlinear effects and MFA in high Mach 
number collisionlcss shocks. 

Ironically, the evidence for large magnetic fields and, 
therefore, nonlinear MFA is obtained exclusively from ra- 
diation emitted by relativistic electrons, while the non- 

1 Physics Department, North Carolina State University, 
Box 8202, Raleigh, NC 27695, U.S.A.; don_ellison@ncsu.edu, 
avladim@ncsu.edu 



linear processes responsible for MFA are driven by the 
efficient acceleration of relativistic ions, mainly protons. 
Relativistic protons can only be directly inferred from 
pion-decay emission at GeV-TeV energies and the un- 
ambiguous interpretation of these observations is difficult 
and remains uncertain. In any case, even if pion-decay 
emission is clearly present, these observations must be 
combined with synchrotron emission from electrons to 
estimate the magne tic field strength. 

A recent report bv lUchivama et all (|2007| ) of rapid time 
variations (~ I yr) in the nonthermal X-ray emission 
from bright filaments in SNR RX JI7I3. 7-3946 provides 
new clues to the electron accelera t ion p rocess in SNRs. 
As emphasized bv lUchivama et all ()2007l ). the rapid time 
variations suggest two conclusions. First, if the varia- 
tions are determined by the synchrotron radiation losses 
of TeV electrons, the magnetic field can indeed be am- 
plified to milligauss-levels at a SNR shock and, second, 
the MFA process is tightly bound with the production 
of relativistic particles (i.e., CRs). An important con- 
sequence of these conclusions is that the fraction of the 
shock ram kinetic energy that is put into accelerated ions 
must be large, e.g., tens of percent. Therefore, DSA must 
be in an efficient, nonlinear regime where the feedback 
of accelerated particles on the shock structure is signif- 
icant. A simple estimate confi r ming this consequence is 
equation (15) in lBell fc Lucekl (|2001h . 

Here we address a single question: Can, as asserted 
by Uchiyama et al. (2007), the large amplified fields 
inferred for electrons from radiation losses in a nonlin- 
ear shock also determine the maximum proton energy 
produced in the SNR shock? We find the answer to 
be no because the inevitable nonlinear shock modifica- 
tion (due to efficient DSA) and the magnetic field vari- 
ation in the shock precursor (due to MFA) make the 
maximum proton energy smaller than what is expected 
without accounting for these effects . Using results from 
IVladimirov. Ellison fc Bvkovi (|2006h for DSA limited by 
the finite size of the shock, we estimate that the max- 
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imum proton energy is at least an order of magnitude 
less than that predicted assuming an unmodified shock 
with a large, nearly constant field inferred from electron 
radiation losses. Our result is similar to that found by 
iBlasi. Amato fc CapriolU (|2007f ) in a time-dependent cal- 
culation of DSA where the acceleration is limited by the 
age of the shock rather than the size, an indication that 
the nonlinear effects we discuss are robust. 

2. ACCELERATION IN A CR MODIFIED SHOCK 

The commonly accepted picture of nonlinear DSA is 
the following: accelerated protons diffuse upstream of 
the shock front, carrying a significant fraction of the 
shock ram kinetic energy, and 'push' on the incoming 
plasma. In order to conserve momentum, the shock forms 
a smooth precursor where the flow speed (in the shock 
frame) gradually decreases from the far upstream speed, 
uq, to a pre-subshock speed, ui, and then sharply de- 
creases in a viscous subshock to the downstream speed, 
U2- Here, i?tot = U0/U2 is the overall compression ratio 
and i? su b = u\/u2 is the subshock compression ratio. If 
MFA accompanies this nonlinear shock smoothing, the 
low ISM field is amplified from Bq = £?ism to B2 on the 
same precursor length scale. 2 

If electrons suffer strong radiation losses, they will 
have a lower maximum energy i<;™ ax than protons and 
will remain in the vicinity of the subshock, spending 
most of their time downstream as they radiate in the 
strongly amplified downstream field i?2- Protons, on the 
contrary, sample the whole precursor and experience a 
weaker mean field (B(x)), which is between -Bism and 
B2 ■ Because the escape of protons from the shock occurs 
far upstream, where the magnetic field is low, one can ex- 
pect that for the most energetic protons, (B(x)) is closer 
to the weak -Bism than to the strong B2, and, conse- 
quently, their maximum energy £?™ ax is lower than would 
be estimated assuming that the shock is unmodified with 
B(x) ~ B 2 everywhere. The latter c ase, B(x) ~ B 2 
everyw here, is implicitly assumed by lUchiyama et al.l 
(2007) when they write t acc = 0.1(S/mG)- 1 ( J B/TeV) yr 
for the proton acceleration time. 

In a size limited shock, the proton maximum energy, 
E™ ax , will be determined when the upstream diffusion 
length of the most energetic protons becomes compara- 
ble to the confinement size of the shock, typically some 
fraction of the shock radius. We model the confinement 
size with a free escape boundary (FEB) at a distance 
Lfeb in front of the shock. Protons that reach this posi- 
tion stream freely away from the shock without produc- 
ing any more magnetic turbulence. Therefore: 

L FEB ~ D(E™*)/u sk , (1) 

where for the diffusion coefficient, D(E), we assume 

D(E) = Ac/3 = V r g c/3 , (2) 

and it s k is the upstream flow speed, A is the scattering 
mean free path, r g is the gyroradius, and 77 is a parame- 
ter that characterizes the scattering strength. None of 
our conclusions depend importantly on 7/ and we as- 
sume 'Bohm' diffusion and set 77 = 1. Since for ultra- 
relativistic protons r g — E/(eB), equations ([1]) and (|2]| 

2 Here and elsewhere the subscript "0" implies far upstream from 
the shock and "2" implies downstream. 



determine -E™ ax for a given Lpgg, i.e., 

£™ aX OC LpEB^sk-Bsk ■ (3) 

For a quasi-parallel, unmodified (UM) shock with no 
MFA, the maximum proton energy -E™ ax |uM is thus de- 
fined by ([3]) with 7i s k-E>sk = uqBq = U0B2, uq being the 
shock speed and B2 being the downstream magnetic field 
derivable from synchrotron emission of accelerated elec- 
trons. However, for a nonlinear (NL) CR modified shock 
of the same physical confinement size, Lfeb, the maxi- 
mum proton energy -E™ ax |NL will be determined by some 
mean value (u(x)B(x)), giving 

^Inl = (u(x)B(x)) 
££ iax |uM u B 2 

For a strongly modified shock, (u(x)B(x)) <C U0B2, 
and in the following we determine (u(x)B(x)) /(wo-B 2 ) 
using the Monte Carlo model de scribed in detail in 
IVladimirov. Ellison & Bvkovl (|2006l ). 

3. MONTE CARLO SIMULATION 

The Monte Carlo model we use (see 
IVladimirov. Ellison fc Bvkovl 120061 . for full details) 
calculates NL DSA and the magnetic turbulence pro- 
duced in a steady-state, plane-parallel shock precursor 
by the CR streaming instability. While the quasi-linear 
approximations we make are only valid in the linear 
regime we nevertheless follow iBell fc Lucekl (|2001| ) and 
assume that the instability doesn't saturate or damp at 
AB w Bq but continues in the nonlinear regime with 
AB ^> B . With this approximation, we are able to 
self-consistently determine the nonlinear shock structure 
[i.e., u(x) vs. x], the MFA [B c g(x) vs. x], and the 
thermal particle injection. 3 

The NL results we investigate do not depend qualita- 
tively on the particular shock parameters as long as the 
sonic Mach number is large enough to result in efficient 
DSA. Here, we use a shock speed = u = 3000 km s , 
sonic Mach number M s m 30, plasma density nigM = 1 
protons cm -3 , and Bq = Bism = 10 /iG, yielding an 
Alfven Mach number, Ma ~ 140. To these parameters 
we add a FEB boundary at Lfeb ~ 0.1 pc, corresponding 
to 10 8 r g! 0j where r ffiU = m v UQc/ {eBo). This size is com- 
parable to the hot spots in SNR RX J1713. 7-3946 and 
produces a proton energy ~ 10 15 eV in our unmodified 
shock approximation. 

Using the above parameters, we simulate two cases: 
a nonlinear solution, where B is amplified from an up- 
stream value Bq = 10 [iG to a downstream value B2 = 
450 /iG (obtained self-consistently by our model), and an 
unmodified solution with a magnetic field set equal ev- 
erywhere to B 2 — 450 [iG. In these two cases we look at 
i?™ ax to see how the prediction of the NL model, conserv- 
ing momentum and energy, compares t o the prediction 
of the UM model, implicitly assumed bv lUchivama et aT] 
(|2007f ). The information about the maximum energy of 
electrons (which are not included in our calculations) can 

3 Note that the Monte Carlo model ignores the dynamic effects 
of electrons and the NL shock structure is determined solely from 
the backpressure of protons and t he magnetic field. While electron 
acceleration can be modeled (e.g., Baring et al. 1999), we only show 
proton spectra here. 
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Fig. 1. — As labeled, the panels show the bulk flow speed in 
units of the shock speed, uq, the effective magnetic field, B e fj, 
[u(x)B c ff (x)\/ (uqB2), and the energy flux, F cn (x), in units of the 
far upstream energy flux, all versus position, x. The subshock is at 
x = and r 9i o = m p uac/ (eBy) (note the change in x-scale between 
linear and log 10 at x = — 10r 9i o)- In all panels, the solid curves are 
the self-consistent, NL result and the dashed curves are for an UM 
shock where energy and momentum are not conserved. For the 
UM shock, the ^-independent B a g = Bo = Bi ~ 450 fiG has been 
chosen to match the amplified field, B2, in the NL result. Note 
that the x-scale in physical units depends on Bo which is different 
in the two cases but that the FEB is set to the same physical 
distance in both cases, i.e., Lpeb ~ 0.1 pc. The drop in energy 
flux at the FEB in the bottom panel results as energetic protons 
escape the NL shock and carry away energy. The total energy flux 
is conserved in the NL case. 

be inferred graphically from the plot of the acceleration 
time (see Fig. [5]). 

Figure Q] shows the shock structure, u(x), the ef- 
fective magnetic field after amplification, B e e(x), and 
u{x)B c s{x)/(uqB2), for the unmodified case (dashed 
lines), and the nonlinear case (solid lines). The bottom 
panel shows the energy flux, normalized to the far up- 
stream value, for the NL case. The smoothing of u{x), 
the weak subshock (i? su b — 2.9), and the increase in i? to t 
above 4 (-Rtot — 9) are clearly present in the top panel for 
the NL case. These three effects must occur to conserve 
momentum and energy if CRs are efficiently accelerated. 
The quantity u(x)B c ff {x)/{uqB2) ~ 0.1 over most of the 
precursor in the NL case. 

In Figure [2] we show the momentum phase-space dis- 
tributions functions, f(p) (multiplied by p 4 ), and the ac- 
celeration time, i acc . The NL effects evident in Fig. [1] 
result in, 

PnT/Pvm < 0.1 , (5) 
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Fig. 2. — The top panel shows proton spectra, p 4 /(p), and the 
bottom panel shows the acceleration time to a given momentum 
for the shocks shown in Fig. [T] In both panels, the solid curve is 
the NL case and the dashed curve is the UM case. The maximum 
momentum in the UM case is about a factor of 10 larger than in the 
NL case. The dotted line in the bottom panel is the synchrotron 
loss time, t syn ~ 12(B/mG)- 2 (-E/TeV) _1 yr, for B = 0.45 mG. 
The NL shock accelerates particles at a slower rate than the UM 
shock. 

and a longer i acc to a given momentum. Here, p^t* 
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jc. We have not at- 



tempted a detailed fit to SNR RX J1713. 7-3946, but 
note that the concave shape of our proton spectrum, 
above the thermal peak , is similar to that obtained by 
iBerezhko fc Volkl (|2006l ) who find a go od fit to the data, 
inclu ding the HESS TeV observations (jAharonian et alJ 
[2006^ . 

The UM result is obtained from the Monte Carlo sim- 
ulation assuming the same "ther mal leakage" model fo r 
injection as in the NL result fe.g.. Uones fc Ellisonlll991f) . 
This injection scheme works self-consistently with mod- 
ifications in the shock structure and overall compression 
ratio, Rtot, to conserve momentum and energy in the NL 
case. In the UM case, the shock structure and i?tot are 
not adjusted and the thermal leakage model produces 
far too many injected particles to conserve momentum 
and energy. For the UM shock to become a test-particle 
shock with energy conservation, far fewer particles would 
need to be injected so that the normalization of the su- 
perthermal f(p) oc p~ 4 power law becomes low enough, 
relative to the thermal peak, so that in contains an in- 
significant fraction of the total shock ram kinetic energy. 
Since we are only interested in comparing p™ ax in the two 
cases, the normalization of the unmodified power law is 
unimportant since Pyivf on ly depends on £feb ■ 

4. DISCUSSION AND CONCLUSIONS 

The possibility of strong MFA in SNR shocks has been 
strengthened by the recent observations of rapid time 
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variability in h o t spo ts in SNR RX J1713.7-3946 by 
Uchivama et al. | (12007). If we accept the conclusions of 
Uchivama et al.l (|2007f ). the ~ 1 yr variations in X-ray 



emission in some hot spots stem from radiation losses 
for electrons and indicate magnetic fields on the order of 
1 mG. Such large fields would almost certainly be caused 
by MFA occurring simultaneously with the efficient pro- 
duction of CR ions in DSA. 

While a number of other interpretations of the X-ray 
and broadband emission in SNR RX J1713. 7-3946 have 
concluded that the magnetic field present in the parti- 
cle acceleration site is considerably less than 1 mG (e.g. , 



siaerably 

Ellison. Slane fe Gaenslerl 120011: lLazendic etHI 12004 
Berezhko fe Volkl l2006t Porter. Moskealenko fe Strone. 
20061) . we have shown that even if the magnetic field 



inferred from e l ectron radiation losses is as high as 
lUchivama et al.l (|2007l ) claim, the underlying physics of 
MFA in DSA shows that this field cannot be simply ap- 
plied to protons to estimate their maximum energy. 

The essential point is that, if MFA to milligauss lev- 
els is occurring as part of DSA, the acceleration must 
be efficient and the system is strongly nonlinear. The 
accelerated particles and the pressure from the ampli- 
fied field feedback on the shock structure (Fig. [1]) and 
this feedback makes the precursor less confining [i.e., 
(u(x)B(x)) <C uqB2\- Therefore, a shock of a given phys- 
ical size will not be able to accelerate protons to an en- 
ergy as large as estimated ignoring NL effects. Even if 
the proton acceleration is limited by the finite age of the 
shock rather that the finite size, as we have assumed, the 
proton £"p nax will be less than in a TP approximation be- 
cause the highest energy protons will have their acceler- 
ation time determined by (u(x) B(x)) rather than 1/0-62- 
This p oint has been made by iBlasi. Amato fe Capriolil 
(2007) using a time-dependent, semi-analytic calculation 
of DSA with MFA. They have derived a general equation 
and determined p^L* f° r a ran ge of Mach numbers and 
choices for the diffusion coefficient when the acceleration 
is limited by the shock age. Considering the differences 
in the models and the specific parameters chosen for our 
examples of young SNRs, the fact that both calculations 
yield PnlV^um ^ 0.1, indicates the robust nature of the 
result which should apply in any shock efficiently produc- 
ing CRs and magnetic turbulence. 

We note that we have neglected the potentially 
important physical effect of magnetic field damp- 
ing in the strong dow nstream turbulence (e.g., 
iPohl. Yan fc Lazar ian 2005) . The importance of d amp- 
ing is still uncertain (see lCassam- Chenai et al]|2007l ). but 
if it occurs it will be most important for high energy pro- 
tons with long diffusion lengths. In effect, damping may 
determine the FEB but shouldn't qualitatively change 
the results for PnlVPum ■ 

Despite the reduction in £™ ax compared to test- 
particle predictions that our results imply, a remnant 
such as SNR RX J1713. 7-3946 might still produce CRs 
up to the knee. The NL example we have presented 
with B 2 — 450 fiG produces protons up to ~ 100 TeV 
in ~ 100 yr in a confinement region of ~ 0.1 pc. If in- 
stead we had taken Lfeb = 1 pc, a size comparable to 
the western shell of SNR RX J1713.7-3946, our NL model 



would produce ~ 1 PeV protons in ~ lOOOyr. Protons 
of this energy are consistent with the ~ 30 TeV 7-rays 
obser ved from SNR RX J1713. 7-3946 (jAharonian et all 
2006) and when the acceleration of heavy ions such as 
Fe +26 is considered, the maximum particle energy ex- 
tends to > 10 16 eV. 

The production of turbulence in efficient DSA 
is an active area of research and the question 
arises if non-resonant wave generation will substan- 
tially modify our results. Fo r ins tance, recently 
iPelletier. Lemoine fc Marcowithl (|2006f ) have investi- 
gated Bell's non- resonant turbulence generation mech- 
anism (|Bel l2004) and found that it can dominate over 
resonant instabilities for high speed shocks (e.g., u s k ~ 
0.1c). They also find, however, that for slower shocks 
of the order we assume he re (m s i< ~ 0.01c), resonant in- 
stabili ties dominate. The IVladimirov. Ellison fc Bvkovl 
(2006|) model assumes resonant instabilities. 

Other effects may be important as well such as parti- 
cles escaping at the FEB. In a size limited shock, the 
escaping energy flux can be large if DSA is efficient. 
For our example, the bottom panel of Fig. [2] shows that 
~ 30% of the energy flux escapes at the upstream FEB 
and this highly anisotropic beam of particles with ener- 
gies near S™ ax moving through the cold ISM should effi- 
ciently generate turbulence. The det ails of this so-called 
"mag netic bootstrap" scenario (e.g., iBlandford fc Funkl 
|2007| ) are just starting to be investigated but such tur- 
bulence might generate higher effective fields with a cor- 
respondingly higher i?™ ax . Despite these possibilites and 
the uncertainties involved in wave generation, we believe 
none of these effects will qualitatively change our con- 
clusions: the precursor field must still range from Bo to 
Bi and the proton maximum energy will be determined 
by a smaller effective field than that estimated from syn- 
chrotron losses for electrons. 

As a final comment we emphas ize a point also made 
bv IBlasi. Amato fc Capriolil ([2007T ). If MFA is occurring 
and the system is highly NL, it may not be possible to ex- 
plain temporal variations in nonthermal X-ray emission 
simply as a radiation loss time. There cannot be varia- 
tions in X-ray emission on short time scales unless the 
accelerator changes in some fashion on these time scales, 
otherwise the radiation would be steady, or varying on 
the shock dynamic timescale, regardless of how short the 
radiation loss time was. Since the injection and acceler- 
ation of protons and electrons is nonlincarly connected 
to the amplified magnetic field, changes in the electron 
particle distribution and changes in the field producing 
the synchrotron emission, will go together and it may be 
difficult to unambiguously determine the field strength 
from temporal variations. 
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